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FeSe with the PbO structure is a key member of the family of new high-T^ iron pnictide and 
chalcogenide superconductors, as while it possesses the basic layered structural motif of edge-sharing 
distorted FeSe4 tetrahedra, it lacks interleaved ion spacers or charge-reservoir layers. We find that 
application of hydrostatic pressure first rapidly increases Tc which attains a broad maximum of 
37 K at ~7 GPa (this is one of the highest Tc ever reported for a binary solid) before decreasing 
to 6 K upon further compression to ~14 GPa. Complementary synchrotron X-ray diffraction at 
16 K was used to measure the low-temperature isothermal compressibility of a-FeSe, revealing an 
extremely soft solid with a bulk modulus, Kq = 30.7(1.1) GPa and strong bonding anisotropy 
between inter- and intra-layer directions that transforms to the more densely packed /3-polymorph 
above ~9 GPa. The non-monotonic Tc{P) behavior of FeSe coincides with drastic anomalies in the 
pressure evolution of the interlayer spacing, pointing to the key role of this structural feature in 
modulating the electronic properties. 

PACS numbers: 74.70.Dd, 74.25.Ha, 61.05. C- 



The a-polymorph of the simple binary FeSe phase has 
recently emerged as a superconductor with an ambient P 
Tc of ~8-13 Its structure comprises stacks of edge- 
sharing FeSe4 tetrahedra with a packing motif essentially 
identical to that of the FeAs layers in the families of 
the FeAs-based high- Tc superconductorsMi^i^i but lack- 
ing any interleaved ion spacers or insulating layers. The 
structural analogy is reinforced by the observation that 
below 70 K the high-temperature crystal structure be- 
comes metrically orthorhombic (space group Cmma)^ 
displaying an identical distortion of the FeSe layers to 
that observed in the iron oxyarsenide familyi^iiS Theo- 
retical calculations also find a very similar 2D electronic 
structure to that of the FeAs- based superconductors 
with cylindrical electron sections at the zone corner and 
cylindrical hole surface sectionsiii Moreover, supercon- 
ductivity in FeSe is very sensitive to defects and dis- 
order and occurs over a limited range of FeSei_5 non- 
stoichiometryji^ 

The effect of applied pressure on Tc provides crucial 
information in differentiating between competing models 
of superconductivity and in the FeSe binary, Tc is ini- 
tially extremely sensitive to P and rises rapidly to 27 K 
at 1.48 GPai^ At the same time, antiferromagnetic spin 
fluctuations present above Tc are strongly enhanced by 
pressure. ^■^ In the FeAs-based superconductors, the re- 
sponse of Tc to pressurization is complex and sensitively 
depends on the composition of the materials and their 
doping level. Both positive and negative initial pressure 
coefficients, dTc/dP have been measured. Typically for 
the REFeAsOi-ajFj; families, dTc/AP is positive at low 



doping levels and switches over to a negative value as 
X increases J^ii^ii^iiiii^ Moreover, for systems where the 
initial dTc/dP is positive, there is a critical value of P 
above which the trend is reversed and Tc then decreases 
upon further pressurization i ^^d^ 

Despite the importance of the evolution of Tc with P 
in understanding the superconducting properties of the 
Fe-based materials, there is currently little information 
on the detailed pressure dependence of their structural 



properties 
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Here we report on the conducting prop- 



erties of FeSe as a function of P up to ^14 GPa; we find 
a very high initial pressure coefficient with Tc maximized 
at 37 K at ^7 GPa. At higher P, dT^/dP becomes neg- 
ative and Tc is reduced to 6 K at ~14 GPa. We also 
study the precise pressure evolution to 12.8 GPa of the 
FeSe structure at 16 K {i.e. well within the supercon- 
ducting regime) by synchrotron X-ray powder diffraction. 
We find that the orthorhombic a-polymorph survives to 
~9 GPa whereupon the structure changes to that of the 
non-superconducting hexagonal /3-phase. Tc initially in- 
creases with P despite the decreasing SeFeSe thickness 
and the increasing distortion of the FeSe4 units while the 
subsequent decrease clearly correlates with changes in the 
FeSe interlayer spacing. 

The FeSe sample used in this work was prepared, as 
reported elsewherCf- Electrical resistivity measurements 
under high pressures were performed by a standard four- 
probe technique using a diamond anvil cell (DAG). Pow- 
dered NaCl was used as pressure-transmitting medium. 
Pressure was applied at room temperature and measured 
by the ruby fluorescence method. The measured values at 
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FIG. 1: (Color online) (a,b) Resistance against temperature 
for FeSe at several applied pressures, (c), (d) Final observed 
(red circles) and calculated (blue solid line) synchrotron X- 
ray (A = 0.41118 A) powder diffraction profiles at 16 K for 
the FeSe sample at (c) 0.25 GPa and (d) 9.0 GPa. The lower 
green solid lines show the difference profiles and the tick marks 
show the reflection positions of the a (upper) and /3 (lower) 
phases, respectively. 



low temperature are used to discuss the pressure depen- 
dence of Tc- The high-pressure synchrotron X-ray diffrac- 
tion experiments at 16 K were performed at beamline 
BLIOXU, SPring-8. The powder sample was loaded in a 
membrane DAC, which was placed inside a closed-cycle 
helium refrigerator. Daphne"^'^ oil was used as a pres- 
sure medium. The applied pressure measured with the 
ruby fluorescence method was increased at 16 K without 
dismounting the cell from the cryostat. The diffraction 
patterns (A = 0.41118 A) were collected using a flat im- 
age plate detector up to 12.8 GPa. Data analysis of the 
diffraction profiles was performed with the GSAS suite 
of Rietveld programs. 

Fig. la shows the temperature dependence of the resis- 
tance of the FeSe sample at various pressures. Supercon- 
ductivity was observed below 13.4 K with a fairly sharp 
onset. Tc is found to increase initially rapidly upon ap- 
plication of pressure with an accompanying increase in 
width and reaches a broad maximum of 37 K at 6.6 GPa 

— for binary solids, this high-Tc is only surpassed by 
CsaCeo (38 K)2^ and MgBs (39 K)3 In a second experi- 
ment, P was increased first directly to 7.6 GPa and then 
to 13.9 GPa where the onset is smaller at 6 K (Fig. 
lb). Pressure release to near ambient and subsequent 
pressurization to 9.7 GPa resulted in smaller values of 

— this may be related to the irreversibility of the a— 
structural transformation (vide infra) (Fig. 2d). 

Inspection of the diffraction profile^^ at 16 K and 0.25 
GPa (Fig. Ic) readily reveals the orthorhombic (O) unit 
cell (space group Cmma) established for a-FeSe below 
70 K at ambient P4 Additional peaks are also evident 
and these can be accounted for by the presence of a mi- 
nority hexagonal NiAs-type /3-FeSe phase (32.4(1)% frac- 
tion). The diffraction datasets collected with increasing 
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FIG. 2: (Color online) Pressure dependence of: (a) the or- 
thorhombic lattice constants a (circles), 6 (diamonds), and c 
(squares), (b) the unit cell volume, V , (c) the [2c/(a-f&)] ra- 
tio, and (d) the superconducting transition temperature, Tc 
of a-FeSe. The solid lines through the data points to 7.5 GPa 
in (a) and (b) are least-squares fits to the Birch-Murnaghan 
equation-of-state, while the dotted lines are guides to the eye. 



P show no structural changes for a-FeSe between 0.25 
and 7.5 GPa with the a//3— phase assemblage remain- 
ing unaltered.— The same structural model was there- 
fore employed in the Rietveld refinements in this pres- 
sure range, revealing a monotonic decrease in the lattice 
constants and unit cell volume with increasing P (Fig. 
2a, b). However, the response of the lattice metrics to P is 
strongly anisotropic with the interlayer spacing showing 
a significantly larger contraction than the intralayer di- 
mensions. As P increases, 2c/{a+b) smoothly decreases 
until at ~4 GPa, it approaches 1. Further increase in P 
leads to an even higher compression of the FeSe inter- 
layer spacing (Fig. 2c) with [2c/{a+b)] = 0.986(2) at 7.5 
GPa. However, as the sample is pressed to 9.0 GPa (Fig. 
Id), the a-phase fraction begins to decrease and then it 
sharply collapses with an almost complete a—>-p transfor- 
mation taking place. At the same time, while the basal 
plane lattice constants continue to contract, c begins to 
increase, resulting in saturation of the volume response 
to pressure and an increasing 2c/{a+b) ratio (Fig. 2a-c). 
Gradual depressurization to ^2 GPa does not lead to re- 
covery of the orthorhombic a polymorph, implying that 
the a— >/3 phase transformation is irreversible at these low 
temperatures. Recovery of a-FeSe (53% fraction) at this 
P necessitated heating of the sample inside the MDAC 
to 300 K,^ 

Fig. 2b shows the pressure evolution of the unit cell 
volume of a-FeSe together with a least-squares fit of its 
16 K equation-of-state (EOS) to the semi-empirical third- 
order Birch-Murnaghan equationi^ The fit results in val- 
ues of the atmospheric pressure isothermal bulk modu- 
lus, Ko= 30.7(1.1) GPa and its pressure derivative, Kq= 
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6.7(6) taking into account the V{P) data to 7.5 GPa. The 
volume compressibility, k — dlnV / dP— 0.033(1) GPa~^ 
implies a soft highly compressible solid. A comparable 
value of the bulk modulus ~33 GPa at 50 K has been 
proposed from neutron diffraction measurements over a 
restricted 0-0.6 GPa range.— The anisotropy in bonding 
of the a-FeSe structure (Fig. 3d) is clearly evident in 
Fig. 2(a) which displays the variation of the orthorhom- 
bic lattice constants with P. a-FeSe is least compressible 
in the basal plane, in which the covalent Fe-Se bonds lie 
(dlna/dP= 0.029(2) GPa'^, dln5/dP= 0.026(3) GPa'^), 
while the interlayer compressibility, dlnc/dP= 0.065(4) 
GPa""'^ is ~2.5 times larger, implying very soft Se-Se in- 
terlayer interactions. 

The immediate consequence of the structural simplic- 
ity of FeSe is that the ~2.91 A-thick SeFeSe building 
blocks are in close proximity separated by each other by 
only ~2.58 A (Fig. 3d)i^ This contrasts sharply with 
the enhanced interlayer separation in the FeAs analogues 
(e.g. interslab separation is ~5.79 A in SmFeAsO which 
comprises interleaved SmO layers^^ and '^3.45 and 3.34 
A in SrFe2As2 and LiFeAs which possess AE^+ and Li+ 
spacers, respectively) As a result, the c axis in FeSe 
is only marginally larger (1.03 times) than the average 
basal plane dimensions. The electronic consequence of 
this structural size proximity should be to render the 
electronic structure of FeSe more 3D in nature than those 
of any of the related FeAs superconductors. Moreover, 
given the softness of the interlayer Se-Se interactions 
relative to the covalently bonded SeFeSe slabs, appli- 
cation of pressure should have a profound influence on 
the structural and electronic dimensionalities, allowing 
their tuning at interlayer contact separations inaccessi- 
ble in other currently known Fe-based superconductors. 
This is of key importance if we recall both the extreme 
sensitivity of Tc to P and its non-monotonic response 
at high P. The distinct structural response to P of a- 
FeSe is clearly apparent by considering the compressibil- 
ity and its anisotropy (Fig. 2a, b). Firstly, the low-T 
bulk modulus, Kq ~ 30.7(1.1) GPa is the smallest mea- 
sured thus far in any of these systems — Kq is 57.3(6) 
GPa for LiFeAs^ 78(2) GPa for LaFeAsOo.gFo.i, and 
102(2) GPa for NdFeAsOo.88Fo.i2'^'^ In addition, the 
compressibility of a-FeSe along the interlayer direction is 
the largest amongst these systems — the c axis contracts 
by 7.3% at 7.5 GPa — while the basal plane compress- 
ibility (contraction by 3.3% at 7.5 GPa) is comparable to 
that in LiFeAs but considerably larger than those of the 
iron oxyarsenides. 

This sets the scene to discuss the pressure response of 
Tc in a-FeSe and its relationship to the structural evolu- 
tion with P. Fig. 3a, 3c show the pressure dependence at 
16 K of the Fe-Se and Se-Se interatomic distances. The 
Fe-Se bonds contract smoothly to ~4 GPa but remain es- 
sentially unchanged with further increase in P resulting 
in an overall decrease of 2.3% at 9.0 GPa. Similarly the 
intralayer Se-Se distances decrease monotonically with a 
somewhat larger contraction of 3.7%. However, the pres- 




FIG. 3: (Color online) (a-c) Pressure dependence of: (a) the 
Fe-Se bond distances, (b) the three distinct Se-Fe-Se bond 
angles, and (c) the inter-layer Se-Se contacts and the height of 
the Se atoms above the Fe plane for the orthorhombic crystal 
structure of a-FeSe at 16 K. (d) Schematic diagram of the 
crystal structure of a-FeSe with selected distances at 16 K 
and 0.25 GPa. (e) Geometry of the FeSe4 units and definition 
of the three Se-Fe-Se bond angles together with their values at 
16 K and 0.25 GPa. (f) Pressure dependence of the [(a-|-6)/2c] 
ratio at 16 K and the superconducting Tc of a-FeSe. The lines 
through the points are guides to the eye. 

sure response of the interlayer Se-Se contacts is much 
steeper with a 9.8% decrease reflecting the very large 
interlayer compressibility — the SeFeSe slabs approach 
each other rapidly up to 7.5 GPa but then their contact 
distance appears to saturate. As a result, there is a cor- 
relation between the pressure evolution of the unit cell 
metrics and the superconducting Tc (Fig. 3f). At low 
pressures, the rapid increase in (a+6)/2c is mirrored by 
the rapidly increasing Tc (a-FeSe has the smallest bulk 
modulus and highest pressure coefficient of Tc amongst 
the Fe-based superconductors). Remarkably the pressure 
range at which the maximum Tc is found coincides with 
the onset of the structural anomalies. While the basal 
plane lattice constants continue to contract smoothly, the 
interlayer spacing begins to expand slightly leading to a 
decreasing {a+b)/2c ratio as P increases above 7.5 GPa. 
Thus it is tempting to ascribe the non-monotonic Tc{P) 
bahavior to the competition between the effects of the in- 
terlayer SeSe interactions (tuning the doping level) and 
the intralayer FeSe bonding and suggests that higher TcS 
are associated with smaller interlayer spacings (increased 
dimensionality) . 

Fig. 3b, c shows the pressure dependence at 16 K of the 
thickness of the SeFeSe slabs and the crystallographic 
bond angles of the FeSe4 tetrahedra (Fig. 3e). It has 
been argued for the iron (oxy) arsenides that the geometry 
of the edge-sharing FeAs4 tetrahedral units sensitively 
controls the width of the electronic conduction band and 
therefore the magnitudes of the slab thickness and the 
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As-Fe-As angles are important parameters in tuning the 
electronic properties of these systems and determining 
^^ 30,31,32 Empirically appears to be maximal when 
the FeAs4 units are close to regular with As-Fe-As angles 
of 109.47°.'^^ Along these lines, the observation of neg- 
ative pressure coefficients of Tc has been rationalized in 
terms of increased tetrahedral distortion away from regu- 
lar shape a-FeSe has an identical distortion mode 
of the FeSe4 tetrahedra with a sizable distortion away 
from regularity (7° angle difference, Fig. 3e). However, 
while the distortion also increases significantly with in- 
creasing P and the SeFeSe slab thickness decreases (Fig. 
3b, c), the pressure coefficient of is positive in this 
range, implying that the empirical generalization which 
appears to hold well for the iron (oxy)arsenides is not ap- 
plicable in the present system — apparently the depen- 
dence of the electronic structure on the interlayer sep- 
aration dominates in determining the superconducting 
properties of a-FcSe. 

Finally, we discuss the high-P transformation. 
Firstly, pressurization leads to a rapid decrease in the 
Se-Se interlayer contacts in the a-phase — at 7.5 GPa, 
these approach a value of 3.3 A, which is significantly 
smaller than the sum of the Se^~ ionic radii. As a re- 
sult, the sign of the compressibilty along c changes and 
the interlayer spacing begins to increase upon further in- 
crease in P. In addition, the three-dimensional /3-FeSe 



is more densely packed than the quasi-two-dimensional 
a polymorph. The large difference in packing density is 
retained at high P and at 9.0 GPa, the a^(3 transition 
is accompanied by a ~14.6% decrease in volume as the 
Fe-Se coordination changes from tetrahedral to the more 
densely packed octahedral one. At the same time, the 
transition results in an increase of the Fe-Fe and Fe-Se 
distances by 6.4% and 4.0%, respectively. These observa- 
tions point toward a scenario where the a-phase cannot 
sustain pressures above 9 GPa as the interlayer separa- 
tion is now so small that is energetically favorable to form 
the more efficiently packed /?-phase with a release of the 
bonding distances. 

In conclusion, we have found that the superconductiv- 
ity onset in a-FeSe attains a broad maximum of 37 K at 
^7 GPa. The orthorhombically distorted superconduct- 
ing phase is extremely soft and the pressure response of 
the structure reveals an intimate link between the SeFeSe 
interlayer separations and the superconducting proper- 
ties. Detailed band structure calculations should be able 
to shed light on the accompanying evolution of the dop- 
ing level of the SeFeSe layers and decipher the relative 
importance between the increased dimensionality and the 
FeSe slab geometry to superconductivity. 
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